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1 Introduction 



There has been sufficient observational indication suggesting a causal connection 
between the binary history of neutron stars and the evolution of their magnetic 
' field. In particular, recent observations and their statistical analyses suggest that 

. (see Bhattacharya 1995, 1996, Bhattacharya & Srinivasan 1995, Ruderman 1995 for 

O ' detailed reviews): 

<N : 

pH . 1. Isolated pulsars with high magnetic fields (~ 10^^ - 10^^ G) do not undergo 

I any significant field decay during their lifetimes; 

' 2. Binary pulsars as well as millisecond and globular cluster pulsars, almost al- 

CN ' ways having a binary history, possess much lower field strengths, down to 
~ 10^ G; 

I 3. Most of these low-field pulsars are also quite old (age ~ 10^ yrs), which implies 

m 



o 



that their field must be stable over such long periods; 



4. Binary pulsars with massive companions, like the Hulse- Taylor pulsar, have 

' field strengths in excess of 10^^ Gauss, whereas the low mass binary pulsars 

O . include both high-field pulsars and very low-field objects like the millisecond 

O ■ pulsars. 

Q^i 5. The evolutionary link between millisecond pulsars and low-mass X-ray bina- 

ries seem to be borne out both from binary evolution models and from the 
comparative study of the kinematics of these two populations. 

^ I Therefore, in this thesis we have tried to understand the mechanism of field evolution 

in neutron stars that are members of binary systems. To this end we have looked 
at four related problems as described below : 

X 

■ • the effect of diamagnetic screening on the final field of a neutron star accreting 

\ material from its binary companion; 

• evolution of magnetic flux located in the crust of an accreting neutron star; 



application of the above-mentioned model to real systems and a comparison 
with observations; 



• an investigation into the consequences of magnetic flux being initially located 
in the core of the star and its observational implications. 



2 Effect of Diamagnetic Screening 

We investigate the effect of diamagnetic screening on the final field of an accreting 
neutron star, in particular we try to answer the following questions : 

1. Are diffusive time-scales, in the layers where the field is expected to be buried, 
long enough to allow screening to be an effective mechanism for long-term field 
reduction ? 
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Figure 1: Variation of the time-scales with surface fields strength. The (a), 
(b), and (c) groups of curves correspond to diffusive, screening and instabil- 
ity time-scales respectively. The curves (in each group) correspond to M = 
10^^*^, 10~^'^, 10~^, 10^^'^ Mq yr~^respectively, the upper curves being of lower val- 
ues of accretion. 

2. Is it at all possible to bury the field or create horizontal components at the 
expense of the vertical one against Rayleigh- Taylor overturn instability ? 

We find that, 

1. The density of flow increases with an increase in the field strength. The larger 
the field, the deeper and denser it gets buried . The flow density also has a 
mild positive dependence on the rate of mass accretion. But even for very large 
values of the surface field strength the flow does not occur at densities beyond 
~ 10^ gcm~^. That means the flow always takes place within the liquid layer. 
And the earlier contentions of a burial within the solid layer does not really 
happen. 

2. As the screening time-scales are always much smaller than the diffusive time- 
scales a condition of flux-freezing prevails and material movement indeed should 
proceed dragging the field lines along. 

3. But, since the instability time-scale (sum of the overturn and the re-connection 
time-scales) is much too smaller than the other two time-scales of the problem, 
(fig. H) any stretching of field line is quickly restored (over the instability time- 
scale) and the material effectively flows past the field lines without causing any 
change. 

Therefore, it is not possible to screen the magnetic field of a neutron star by the 
accreting material in order to reduce the magnitude of the external dipole observed. 

3 Evolution of Crustal Magnetic Field in an Accreting 
Neutron Star 

A possible mechanism of field decay is that of rapid ohmic diffusion in the accretion 
heated crust. The effect of accretion on purely crustal fields, for which the current 
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Figure 2: Evolution of the surface magnetic field for six values of accretion rate. 
The curves 1 to 5 correspond to M = 10"^^ lO'^^ 10-i\ 10"^°, 2.0 x 10-^° Mq 
yr~^with the crustal temperatures obtained from logT = 0. 397 log M + 12.35. The 
dashed curves 6a, 6b and 6c correspond to T = 10^'^, 10^'^^, 10^'^ K respectively 
for an accretion rate of M = 10~^ Mq yr~^. All curves correspond to Q (impurity 
content of the crust) = 0.0, but are insensitive to the value of Q. 

loops are completely confined within the solid crust, is two-fold. On one hand the 
heating reduces the electrical conductivity and consequently the ohmic decay time- 
scale inducing a faster decay of the field. At the same time the material movement, 
caused by the deposition of matter on top of the crust, pushes the original current 
carrying layers into deeper and denser regions where the higher conductivity slows 
the decay down. The mass of the crust of a neutron star changes very little with 
a change in the total mass; accretion therefore implies assimilation of the original 
crust into the superconducting core. When the original current carrying regions 
undergo such assimilation, further decay is stopped altogether. 

The combination of enhanced ohmic diffusion due to crustal heating and the trans- 
port of current-carrying layers to higher densities due to the accreted overburden, 
causes the surface field strength to exhibit the following behaviour (fig. ^): 

1. An initial rapid decay (power law behaviour followed by exponential behaviour) 
followed by a leveling off (freezing), 

2. Faster onset of freezing at higher crustal temperatures and at a lower final 
value of the surface field, 

3. Lower final fields for lower rates of accretion for the same net amount of ac- 
cretion, 

4. The longer the duration of the pre-accretion phase the less the amount of field 
decay during the accretion phase, and 

5. The deeper the initial current loops are the higher the final surface field. 
For details see Konar &: Bhattacharya (1997). 
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Figure 3: Evolution of the surface magnetic field in LMXBs with an wind accretion 
rate of M = 10"^'* M0 yr~^. Curves 1 to 5 correspond to initial current config- 
urations centred at p = 10^^, 10^^'^, 10^^, 10^^'^, lO^^g cm~^. All curves correspond 
to Q = 0.0. A standard cooling has been assumed for the isolated phase. The 
wind and Roche-contact phases are plotted in expanded scales. The dashed and 
the solid curves correspond to accretion rates of M = 10~^, 10"^*^ yr~^in the 
Roche-contact phase. 

4 Comparison with Observation 

We have investigated the evolution of the magnetic field of neutron stars in its en- 
tirety - in case of the isolated pulsars as well as in different kinds of binary systems, 
assuming the field to be originally confined in the crust. We model the full evolution 
of a neutron star in a binary system through several stages of interaction. Initially 
there is no interaction between the components of the binary and the evolution of 
the neutron star is similar to that of an isolated one. It then interacts with the 
stellar wind of the companion and finally a phase of heavy mass transfer ensues 
through Roche-lobe overflow. 

We find that the model can explain almost all the features that have been observed 
to date. And our conclusions can be summarized as follows : 

• for this model to be consistent with the statistical analyses performed on the 
isolated pulsars at the most a maximum value of 0.05 for the impurity strength 
can be allowed; 

• HMXBs produce high-field long period pulsars provided the duration of the 
wind accretion phase is short or the initial current distribution is located at 
higher densities; 

• Relatively low-field {B 10^" Gauss) objects near death-line (low-luminosity 
pulsars) are also predicted from HMXBs; 

• LMXBs will produce both high-field long period pulsars as well as low-field 
short period pulsars inclusive of millisecond pulsars in the later variety (fig. 
^; and 

• a positive correlation between the rate of accretion and the final field strength 
is indicated that is supported by observational evidence. 
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Figure 4: Evolution of the surface magnetic field in low mass X-ray binaries. The 

dotted and the solid curves correspond to accr etion rates of M = 10"^ 10"^° Mq 

yr~^in the Roche contact phase. The curves 1 to 5 correspond to Q = 0.0, 0.01, 

0.02, 0.03 and 0.04 respectively. All curves correspond to a wind accretion rate of 
M = 10-16 yj.-i_ 

For details see Konar & Bhattacharya (1999a). 

5 Spin-down induced Flux-expulsion and its Consequences 

Finally, we look at the outcome of spindown-induced expulsion of magnetic flux 
originally confined to the core, in which case the expelled flux undergoes ohmic decay. 
The general nature of field evolution again seems to fit the overall observational 
scenario. The nature of field evolution is quite similar to that in case of a purely 
crustal model of field evolution though the details differ. Most significantly, this 
model has the requirement of large values of the impurity strength Q in direct 
contrast to the crustal model. To summarize then : 

• The field in isolated neutron stars do not undergo any significant decay, over 
the active lifetime of the pulsar, conforming with the statistical analyses. 

• The field values in the high mass X-ray binaries can reaming fairly large for a 
moderate range of impurity strength. 

• A reduction of three to four orders of magnitude in the field strength can be 
achieved in the low mass X-ray binaries provided the impurity strength is as 
large as 0.5 (fig. ^). 

• If the wind accretion phase is absent then to achieve millisecond pulsar field 
values, impurity strength in excess of unity is required. 

For details see Konar & Bhattacharya (1999b). 

6 Conclusions 

In this thesis we have mainly investigated two models of field evolution - that of 
an initial crustal current supporting the field and that of spin-down induced flux 
expulsion, also looked at the effect of diamagnetic screening in an accreting neutron 
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star. Following are our conclusions regarding the nature of field evolution in the 
crust. 

• Pure Ohmic Decay in Isolated Neutron Stars : 

1. A slow cooling of the star gives rise to a fast decay and consequent low 
final field. The opposite happens in case of an accelerated cooling. 

2. An initial crustal current distribution concentrated at lower densities 
again gives rise to faster decay and low final surface field. Whereas if 
the current is located at higher densities the decay is slow resulting in a 
higher final surface field. 

3. A large value of impurity strength implies a rapid decay and low final 
field. If the crust behaves more like a pure crystal the decay slows down 
considerably. 

• Accretion-Induced Field Decay in Accreting Neutron Stars : 

1. In an accreting neutron star the field undergoes an initial rapid decay, 
followed by slow down and an eventual freezing. 

2. A positive correlation between the rate of accretion and the final field 
strength is observed, giving rise to higher final saturation field strengths 
for higher rates of accretion. 

3. An expected screening of the surface field by the diamagnetic accreting 
material is rendered ineffective by the interchange instabilities in the liq- 
uid surface layers of the star. 

4. To produce millisecond pulsars in LMXBs in spin-down induced flux ex- 
pulsion model very large values of impurities are required. This makes the 

surface field go down to very low values in 10^ years in isolated pulsars 
in contrast to a purely crustal model. 
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